Quantum chaos and fluctuations in isolated nuclear spin systems 
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Using numerical simulations we investigate dynamical quantum chaos in isolated nuclear spin 
systems. We determine the structure of quantum states, investigate the validity of the Curie law 
for magnetic susceptibility and find the spectrum of magnetic noise. The spectrum is the same 
for positive and negative temperatures. The study is motivated by recent interest in condensed- 
matter experiments for searches of fundamental parity- and time-reversal-invariance violations. In 
these experiments nuclear spins are cooled down to microkelvin temperatures and are completely 
decoupled from their surroundings. A limitation on statistical sensitivity of the experiments arises 
from the magnetic noise. 

PACS numbers: 05.50.-|-q, 05.45.Pq 
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Introduction. Much of the present knowledge about 
violations of the fundamental symmetries P (invariance 
under spatial inversion) and T (invariance with respect to 
time reversal) comes from experiments measuring P- and 
T-violating permanent electric dipole moments (EDM's) 
of atoms, molecules and the neutron, see, for example, 
Ref. 0. Most EDM experiments measure precession 
of the angular momentum of the system in an applied 
electric field analogous to the Larmor precession in an 
applied magnetic field. 

In addition to such precession experiments, there are 
EDM searches of another kind 2,^]jwhich have drawn 
recent renewed attention 0, H S i3 fi II E3 ■ The idea 
of these experiments is the following. Suppose that we 
have some condensed matter sample with uncompensated 
spins. If an electric field is applied to the sample, it in- 
teracts with the associated (P- and T-violating) EDM's, 
leading to a slight orientation of the spins in the di- 
rection of the electric field. This orientation, in turn, 
is measured by measuring the induced magnetization of 
the sample. In this work we concentrate on effects re- 
lated to nuclear spins in insulators with fully compen- 
sated electron spins 0, 1^ ■ The EDM energy shift under 
discussion is about 10~^^ — 10~^^eV per nuclear spin. 
For comparison a similar shift is created by a magnetic 
field B - 10-16 - 10-20T. So the effect is tiny and hm- 
itation to statistical sensitivity comes from the fact that 
the number of spins, n ~ 10^^, inspite of being large is 
still finite Q . Basically the limitation comes from a kind 
of magnetic shot noise. In the EDM experiments the nu- 
clear spins must be cooled down at least to 100/iK and 
optimally down to 10— lOOnK. At low temperatures the 
spins are completely decoupled from the crystal lattice, so 
there is no contact with any heat bath. This motivates 
the problem considered in the present work: magnetic 
noise of an isolated nuclear spin system. An important 
point is that the total spin of all nuclei is not conserved 
because magnetic dipole-dipole interaction depends on 
relative orientation of nuclei. 

Concerning previous work we first of all refer to the 
1959 paper by Hebel and Slichter where the physical 



meaning of temperature for an isolated spin system has 
been discussed. This work assumed the validity of a sta- 
tistical approach for the isolated quantum system. This 
is what nowadays is called dynamical quantum chaos. 
The problem of the onset of dynamical quantum chaos 
has been addressed much later. Level statistics in a spin 
system with a mobile fermion has been investigated by 
Montambaux et al in 1993 The criterion for on- 

set of quantum chaos in spin glass shards with Heisen- 
berg interaction in a random external magnetic field has 
been derived by Georgeot and Shepelyansky in 1998 0|. 
There have been also numerical studies of level statistics 
in XYZ spin chains with and without a random magnetic 
field that demonstrated non-Wigner behavior in the ab- 
sence of magnetic field Q . 

Model. To be specific, we consider the Lead Titanate 
ferroelectric suggested for EDM experiments in Refs. 
0, ■ The interaction between the nuclei is governed 
by the magnetic dipole-dipole interaction, the strength 
of which falls off as 1/r^, Jap = 7^(^a^ — 'inaUpj/r^ . 
Therefore in the paramagnetic phase that we are inter- 
ested in it is sufficient to only include nearest-neighbour 
interactions. The dominating contribution to the dipole- 
dipole interaction is from the ^°^Pb isotope which has 
a natural abundance of 22.1% and are distributed ran- 
domly through the lattice. As a consequence the inter- 
action is anisotropic and random in strength. Therefore 
we adopt a model Hamiltonian of the form. 
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where J^Jj is the interaction between spins S = 1/2 on 
nearest sites k and / on a lattice and B is a uniform ex- 
ternal magnetic field. The interaction J^'^ is represented 
by random numbers uniformly distributed between [-J,J]. 
The typical value of J is J ~ lO'^'^eV ~ lOnK 0. For 
the sake of simplicity unless otherwise is stated we take 
a usual square lattice with periodic boundary conditions. 
We will consider the case when the tensor Jaf3 has only 
diagonal components, a — (3, and the case when it has 
both diagonal and off diagonal components. The latter 
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case corresponds to real dipole-dipole interaction. We 
diagonalize the Hamiltonian |^ numerically exactly and 
we need to know all eigenstates and all eigenenergies, so 
the size of the matrix is 2" where n is the number of spins. 
Therefore, practically we are able to consider only rela- 
tively small clusters with n=8,10,12. For a sufficiently 
large lattice all results are expected to be self-averaged. 
However, we consider relatively small clusters, therefore, 
to improve statistics we average results over 100 random 
realizations of J^^. We do not consider an odd number of 
spins to avoid Kramers degeneracy of spectra that cannot 
be important in the thermodynamic limit. 

Structure of chaotic eigenstates. Small magnetic field. 
At zero magnetic field the Hamiltonian l(T]l is invariant 
under time reversal (T-reversal) . This implies that the 
eigenfunctions of the Hamiltonian are split in two sectors 
of different time-parity: in the first sector they stay the 
same and in the second one they change sign (-) un- 
der the action of all the spins being flipped. Therefore 
the expectation value of magnetization that is a T-odd 
operator vanishes for any state, {i\Sz\i) = 0. Because 
of random J, chaos is established in each of these two 
sectors but the sectors do not interact. Anderson local- 
ization of single particle states in a random potential is 
a well known effect. In principle a spatial localization 
of many-body quantum states is also possible, while it 
is known that interaction tends to destroy Anderson lo- 
calization [l^ . Spatial localization would imply Poisson 
level statistics (the distribution of level spacing between 
closest levels), Pp{s) — exp(— s), within a sector with 
a given T-parity in a sufficiently large system. On the 
other hand the Wigner-Dyson distribution Pwd{s) — 
^ exp(— 7rs^/4) within a given sector indicates a full 
chaotization including delocalization. To characterise 
to what degree the statistics of levels are Wigner 
or Poisson, following Ref. we use the parameter 

V = Io° [m - PwDis)] ds/ j;° [Ppis) - PwDis)] ds , 

where P{s) are the statistics measured in numerical simu- 
lations and So = 0.4729 is the intersection point of Pp{s) 
and PwD- Localization effects are always enhanced in 
the ID case. Therefore to investigate the localization sce- 
nario we studied level statistics within a given T-parity 
sector for the Hamiltonian iQJ on a ID ring at -B — 0. 
We found that r/ = 0.18,0.051,0.035 for n=8, 10 and 12 
respectively. So, we do not observe any deviation from 
the Wigner-Dyson distribution within the accessable sys- 
tem size. Thus, we come to a conclusion that there is no 
spatial localization of quantum states in the random spin 
system. This conclusion is similar to that of Ref. 0| for 
mobile interacting particles. 

The combined statistics of levels including both sectors 
at 5 = is given by the sum of two Wigner distributions, 
giving an intermediate statistics with rj « 0.5, see FigQ 
A very small critical magnetic field is needed to mix the 
sectors and hence to lead to a single Wigner distribution. 
The critical magnetic field Bd is given by the condition 



that mixing of two nearest states from different sectors 
is of the order of unity 
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Here 



and are opposite T-parity eigenstates 



of the Hamiltonian at zero magnetic field. The 
level spacing is roughly equal to AE ~ J-y/n/2"', where 
the factor y/n comes from the total width of the spec- 
trum that is discussed below. To estimate a typi- 
cal mixing matrix element in ((21 we use the sum rule 
E,{^+\S.\j-){j-\S.\i+) = {^+\S'A^+) - {i+\EmS'.m + 
12m^k Sz,mSz,k\i+) ~ J ■ Here we take into account that 
J2m ^Im = '^/^ and that X^m^^fc ~* ^ because it is inco- 
herent. There are 2""^ terms in Therefore the typ- 
ical mixing matrix element is |(j-|>5'2|i+)| ^ \/ n/2"+^. 
Hence substitution in (|2Jl gives the following estimate 
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To confirm this analytical estimate we present in Fig. ^ 
the value of rj calculated numerically for different mag- 
netic fields. As we already mentioned at zero magnetic 
field r] « 0.5 indicating intermediate statistics, and at 
large field 77 ^ indicating the Wigner-Dyson distribu- 
tion. Taking the value rj = 0.25 as a crossover point we 
find Bc{n + 2)/Bc{n) « 0.5 in good agreement with Eq. 
©• Thus the value of the first critical field Bd drops 




FIG. 1: r; as a function of the magnetic field: solid line n ■ 
dashed line n — 10, dotted line n = 12. 



down exponentially with number of spins. For the EDM 
experiment we are interested in large systems, n ~ 10^^. 
So in this situation B^i = 0. Note that Eq. Q is not a 
criterion for onset of quantum chaos. This is the crite- 
rion for a crossover from a two component chaotic dis- 
tribution to a single component chaotic distribution or 
in other words the criterion for destruction of T-parity. 
This is quite similar to the dynamical enhancement of 
spatial parity violation in nuclei . 

"Strong" magnetic field. In a strong uniform mag- 
netic field stationary states can be classified by total 
spin projection on the direction of the field, Sz- This 
regime is realized when the level spacing B is larger 
than the matrix element between direct ly c oupled states 
that is of the order of J, see e. g. [iSj- Thus the 
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crossover to the strong magnetic field regime happens 
when B > Bc2 ^ J 10~^T (the numerical value cor- 
responds to Lead Titanate). In the strong field regime 
due to random exchange interaction J the states with a 
given value of Sz are completely mixed up and their ener- 
gies, according to our calculations, are spread within the 
band of width O.Sy^J. Note that "conservation" of 
total Sz does not mean that the z-projection of a partic- 
ular spin is conserved. Every particular spin fluctuates 
very strongly. Thus the energy spectrum consists of suc- 
cessive bands, separated by B, the width of the band is 
about ~ O.Sy^J. There is the Wigner-Dyson statistics 
of levels within a given band. When n » (B/J)'^ the 
bands overlap and this leads to Poisson level statistics. 

The above discussion of the structure of quantum 
states is valid in the case when the tensor J^^ in Hamilto- 
nian JQ) has both diagonal and off diagonal components. 
In the case without off diagonal components there is a 
hidden integral of motion. The Hamiltonian acting on a 
basis state with a particular Sz only changes Sz by 0, ±2, 
therefore Sz{mod2) is conserved. This means that eigen- 
functions can be expanded in a basis with either even 
or odd Sz components. Then the above considerations 
are valid separately for 5'z-even and for Sz-odd sectors. 
It is possible that the non-Wigner behaviour observed in 
Ref. [lj| is due to the hidden integral of motion. 

Temperature, average energy and magnetic susceptibil- 
ity. The textbook analysis of a spin system based on as- 
suming the validity of the Boltzmann distribution gives 
the following well known relations 

X=^, E{T)^E,-^(Ef-Ei) . (4) 

Here x is the magnetic susceptibility, E{T) is the aver- 
age energy at a given temperature (we set the Boltzmann 
constant equal to unity), Ei is the average energy of sta- 
tionary states and Ef is the average energy squared. Ei 
and Ef are independent of temperature. First we want 
to check the validity of Eq.Q) for the isolated dynam- 
ical system In this case Eq.Q in essence defines 
an effective temperature, see the discussion in Ref.|ll(. 
To give a precise meaning to Eq. Q one has to con- 
sider many energy levels inside a bin around some given 
energy E. Then, according to Q the susceptibility av- 
eraged over these levels is related to the energy E. We 
have checked numerically that for the Hamiltonian 
Ei and E^ w O.lSJ^n. Hence, according to 
X ~ — 1.94£'/J^. This agrees well with results of numer- 
ical simulations with Hamiltonian ^ shown in Fig. 
This figure represents values of the level magnetization 
{i\Sz\i) calculated in the field B = O.IJ, averaged over 
bins of width AE = 0.5J, and averaged over random 
Jai3- Deviation from linear dependence at large positive 
and negative energy is the finite size effect: Eqs. (@J) make 
sense only if the average energy is smaller than the "band 



width", \E\ < y/nJ. It is instructive to consider instead 



3.5 




-4 -2 2 °4 '■"''O 0.2 0.4 0.6 0.8 1 

E/J B/J 
FIG. 2: a. Magnetization {i\Sz\i) averaged over energy bins 
AE = 0.5 J. The magnetic field B = O.IJ. Solid, dashed, 
and dotted lines correspond to clusters of size n — 8, 10, 12 
respectively. 

b. The value TM/B versus magnetic field for temperatures 
T=-3J (solid), -J (dashed), J (chain), 3 J (dotted) and for 
cluster sizes n = 8, 12. The lower group of lines correspond 
to n = 8 and the upper group to n = 12. 

of "rectangular binning" another way of sampling of en- 
ergy levels. In particular the Boltzmann sampling defined 
as M = i E^{^\S^\^)e-^'/^, where Z = J2^ e"^'/^. Let 
us stress that temperature here is just a parameter that 
samples a set of energy levels around the average en- 
ergy that according to 10} is = — 0.13nJ^/T. One can 
also check that the energy "window" around the aver- 
age energy is SErms = vOT3nJ. Note that the Boltz- 
mann sampling does not imply a validity of the equi- 
librium Boltzmann distribution. For example we can use 
Boltzmann sampling in a system without interaction that 
never comes to an equilibrium. The value of TM/B cal- 
culated for different cluster sizes, temperatures (positive 
and negative) and magnetic fields and averaged over ran- 
dom JajS is presented in Fig. ^p. For high temperatures, 
|T| = 3J, the results agree perfectly with the Curie law. 
Unexceptedly it agrees even at -B < Bd- For lower tem- 
peratures, \T\ ~ J the results deviate from Curie's law as 
one should expect. To complete the analysis we have also 
calculated the dispersion of the susceptibility and found 
that Sx = Vix^) - (X)' - O.OIV^/T. 

Magnetic susceptibility at nonzero frequency, fluctua- 
tions. There are no temporal fluctuations in an isolated 
system in a given quantum state. The fluctuations come 
from the fact that a large system cannot be prepared in 
a flxed quantum state. We always deal with a density 
matrix that represents a combination of a large number 
of quantum states around some average energy. In this 
section we consider only the Boltzmann sampling: the 
average energy is .E = — O.lSrtJ^/T and the energy "win- 
dow" is SErms — VO.lSrtJ. In this case we can use the 
standard formula for the imaginary part of the magnetic 
susceptibility 0] 

X \{t\Sz\j)\^6{uj + uj,j) . (5) 
Here — (Ei — Ej)/h. First we perform direct numer- 
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ical simulations using this formula and quantum states 
generated by the Hamiltonian ^ at B = O.IJ 0| 
For simplicity we consider the case when the tensor Jq,^ 
has only diagonal components. Results of simulations 
for cluster n = 12 are shown in Fig. 13 The plots for 
71 = 8, 10 are very similar. The result unambiguously in- 
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FIG. 3; TImx/n versus u. Solid line, analytic. Numerical 
results for n=12atT = 3 and T — 1 are shown by dashed 
and dotted lines respectively. 




noise is independent of temperature and moreover it is 
the same for positive and negative temperature. What is 
usually called the fluctuation dissipation theorem in the 
case of negative temperature becomes the fluctuation ra- 
diation theorem because Im[x{uj)] changes sign. 

We have investigated the structure of chaotic quantum 
states in a spin lattice system with random interactions. 
We also checked the validity of the Curie law for magnetic 
susceptibility and find the spectrum of magnetic noise. 
The temperature independent noise limits the statistical 
sensitivity of experiments on parity and time-reversal- 
invariance violations. 

We thank D. Budker. J. Imry, A. Luscher, J. Oitmaa, 
and A. Sushkov, for stimulating discussions and critical 
comments. This work was supported by the Australian 
Research Council. 



dicates a Gaussian line shape, Im[x] « ■^y/TTUJTe~'^ , 
with r « 0.85/ J. Assuming the Gaussian shape we 
can also calculate r analytically using a method usually 
applied in NMR studies At a large temperature, 

Ei <^ T, and hence we can represent Q as Im[x{uj)] = 
^f{Lj) where the lineshape /(w) is given by, /(w) = 
J2i j l(*l'5'zb)P'^(^ + ^ij)- Let us calculate the second 
moment of /(w), M2 = J^^ uj'^ f {uj)duj ^ /(w)dw. 
Using the completeness relation and the spin decoupling 
= J2ki^^kSzi ^/^J2ki^ki (the indexes k, I enu- 
merate sites) we find the denominator /(w)dw = 
n2"~^. To calculate the numerator one represents it in 
terms of the commutator of spin with the Hamiltonian 
/.!l^V(^^)dw = J2M[H,Sz]^\i). Then calculating the 
commutator and using the spin decoupling wc find 



/Hd^ = -^[(j2^ + j2^ + j2j 



(6) 



(ki) 



{JxxJyy + JyyJzz + JxxJzz) + '^{Jxy + ^zy + ^xz)] ■ 

Here {kl) represents a pair of nearest sites. Since we 
have performed numerical simulations without off diago- 
nal components of J^^, the right hand side of Eq. © is 
^2"-i J2/3 and hence M2 = 2/3 J^. On the other hand 
for the Gaussian line shape M2 = l/(2r^). Hence the 
analytical calculation gives the r = ■\/3/4/ J. The Gaus- 
sian curve with this value of the relaxation time is shown 
in Fig. El by solid line. It is in good agreement with re- 
sults of direct numerical simulations. The real part of the 
susceptibility can be easily found using Kramers-Kronig 
relations (see, for example, Ref. |[18j|) 

As soon as is known then using the fluctuation 

dissipation theorem one can find magnetic noise, i. e. 
the spectral density (M^);^ of the square of the devi- 
ation of the magnetization from its average value 0, 
V^{M^)^ = hcoth{huj/2T)Im[x{uj)]. V is volume of the 
sample. A very interesting point is that at w <C T the 
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